Atom-dimer scattering and long-lived trimers in fermionic mixtures 
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We consider a heteronuclear fermionic mixture on the molecular side of an interspecies Feshbach 
resonance and discuss atom-dimer scattering properties in uniform space and in the presence of 
an external confining potential, restricting the system to a quasi-2D geometry. We find that there 
is a peculiar atom-dimer p-wave resonance which can be tuned by changing the frequency of the 
confinement. Our results have implications for the ongoing experiments on Lithium-Potassium 
mixtures, where this mechanism allows for switching the p-wave interaction between a K atom and 
Li-K dimer from attractive to repulsive, and forming a weakly bound trimer with unit angular 
momentum. We show that such trimers are long-lived and the atom-dimer resonance does not 
enhance inelastic relaxation in the mixture, making it an outstanding candidate for studies of p- 
wave resonance effects in a many-body system. 



PACS numbers: 

The last several years have seen enormous advances in 
the field of ultracold fermionic gases. The main progress 
has been made in studies of the crossover from the 
Bardeen-Cooper-Schrieffer superfluidity to Bose-Einstein 
condensation of molecules (BCS-BEC crossover) in mix- 
tures of two hyperfine states of the same atom (either Li 
or K) close to an intercomponent Feshbach resonance pQ . 
Tunability of these systems has allowed one to investigate 
the effect of a population imbalance on the BCS-BEC 
crossover [2J [3] , study vortices in strongly interacting ro- 
tating fermionic supcrfluids 4 , and explore other exotic 
phenomena. The BCS-BEC crossover in a mixture of 
different fermionic species and the formation of weakly 
bound heteronuclear bosonic dimers is now the subject 
of very active investigation. Naturally, most of the ex- 
perimental progress has been made for the alkali mixture 
of 6 Li and K. In particular, Taglieber and co-workers 
have cooled the system to quantum degeneracy [5] , inter- 
species Feshbach resonances have been studied in a non- 
degenerate system [Bj |7] , and, recently, the formation of 
weakly bound Li-K dimers at two of these resonances 
has been reported [5J [5] . The association of atoms into 
weakly bound dimers is an important step toward ob- 
taining ultracold dipolar molecules, as has recently been 
demonstrated for fermionic 87 Rb- 40 K molecules at JILA 
[10] . The same technique can be used to transfer bosonic 
6 Li- 40 K dimers to the ground rovibrational state with a 
large dipole moment |llj . 

Although the dipole moment of weakly bound hete- 
ronuclear dimers is negligible, they can strongly interact 
by exchanging light atoms. For sufficiently large mass ra- 
tios this exchange interaction dominates over the kinetic 
energy of the dimers and can even lead to their crystal- 
lization [T^] • On the few-body level the role of the mass 
imbalance is distinctly visible considering the problem of 
two identical heavy fermions of mass M interacting res- 
onantly with a light atom of mass to. For mass ratios 



M/m > 13.6 an attractive effective potential mediated 
by the exchange of the light atom dominates over the 
centrifugal barrier for the identical heavy fermions. This 
leads to the Efimov effect - the existence of an infinite 
number of bound trimer states [13) . For smaller mass 
ratios the centrifugal barrier is dominant. It excludes 
the Efimov effect and also suppresses recombination pro- 
cesses requiring three atoms to approach each other to 
very short distances, which is crucial for the collisional 
stability of the gas |14j . Amazingly, M/m — 13.6 is not 
the only critical mass ratio for this system. Recently, 
Kartavtsev and Malykh [15, have shown that for positive 
interspecies scattering length a and for M/m > 8.2 there 
is a weakly bound not Efimovian trimer state, which for 
smaller mass ratios turns into a p-wave atom-dimer scat- 
tering resonance. 

In this Letter we consider a heteronuclear fermionic 
mixture with M/m < 8.2 and study the scattering of a 
heavy atom by a weakly bound heteronuclear molecule, 
pointing out not only the practical relevance but also 
uniqueness of the 6 Li- 40 K case (M/m as 6.64). We cal- 
culate s- and p-wave atom-dimer scattering phase shifts 
as functions of the collision energy taking into account 
the finite width of the interspecies Feshbach resonance. 
We show that for sufficiently small dctunings the p-wave 
atom-dimer interaction is resonant and should dominate 
the dynamics of the atom-dimer mixture in the ultra-cold 
regime. We also show that the atom-dimer p-wave reso- 
nance can be tuned by confining the system to a quasi- 
2D geometry, and one can even turn it into a real trimer 
state. Finally, we predict that the resonant p-wave atom- 
dimer interaction does not amplify the collisional relax- 
ation in the system. In particular, the lifetime of a trimer 
can exceed several seconds. The relaxation in s-wave 
atom-dimer and dimer-dimer collisions, more dangerous 
in the case of a narrow Feshbach resonance, can also be 
suppressed by decreasing the magnetic field detuning. 
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In order to solve the three-atom problem at hand we 
use the zero-range model, which is well justified for ul- 
tracold gases where the van der Waals range of the in- 
teratomic potential, R e , is orders of magnitude smaller 
than all other lengthscales. The method has been ap- 
plied to three unequal-mass fermions in free space near 
a wide Feshbach resonance [TS]. The effect of the finite 
width of the resonance is described in Ref. [17] in the 
case of three identical bosons. Mora and co-workers [TS] 
have included an external quasi- ID confinement in the 
treatment of the three equal-mass fermion problem. The 
synthesis of these techniques is outlined below. 

Consider a pair of identical heavy fermions of mass 
M interacting with a light atom of mass m via an s- 
wave Feshbach resonance. We assume that the atoms 
are harmonically confined in the z-direction with the sin- 
gle particle Hamiltonians Hi = -V^./2M + MQ 2 Z 2 /2 

and h = —W 2 /2m + mu} 2 z 2 /2, where we set h = 1. 
In the zero-range approximation the interaction acts 
only on the "boundaries" r = Ri and r = R2 leading 
to the l/|r — Ri ^singularities in the 3-body wavefunc- 
tion ^(Ri, R2, r). We introduce an auxiliary function 
/(Ri,R 2 ) = lim r ^ Rl 47r|r-Ri|*(Ri,R 2 ,r), which, in 
fact, contains all the information about the state of the 
system as \& can be retrieved by solving the equation 



(ifi + H 2 + h - E)V 



/(Ri,R 2 )(S(r-Ri)/2^ 
-/(R 2) R 1 )5(r-R 2 )/2 M , (1) 



where /1 = mM/(m + M) is the reduced mass. Below the 
three-atom threshold, E < D, + u>/2, the mapping / — > "J 
is unambiguous and is given by an integral operator, the 
kernel of which is obtained from the Green function of 
the left hand side of Eq. 0. In the limit r — > Ri this 
integral expression for ^ reduces to the form 

*(Rx,R 2 ,r) - (l/|r - R x | - £)/(R x , R 2 )/47r, (2) 

where the integral operator K gives a regular contribu- 
tion in the considered limit. 

The singular and regular terms in Eq. ^ should 
match the two-body scattering problem and, hence, are 
related by the Bethe-Peierls boundary condition W oc 
(l/|r — Ri| — l/a(E c )), where the energy dependent scat- 
tering length defined as a(Ec)^ 1 = a^ 1 + 2p,R*E c is eval- 
uated at the collision energy of the two atoms [T7] . Here 
a is the scattering length (at zero collision energy) and 
the length parameter R* is inversely proportional to the 
width of the Feshbach resonance. The collision energy, 
E c , can be obtained by separating the relative motion of 
the colliding pair and subtracting the energies of all other 
degrees of freedom from the total energy of the system. 
Thus, E c can be written as a differential operator acting 
on /: 



E e f = 



E 



(Mft 2 + mto 2 )Z 2 



2(M- 



H 2 



Finally, the integro-differential equation for / reads 

[-K + a- 1 +2/iiT£y/(Ri,R 2 ) = 0. (3) 

Equation ^ can be regarded as a two-body Schrodinger 
equation describing the motion of a dimer and an atom, 
/(Ri,R 2 ) being the corresponding wavefunction. 

Let us now consider the atom-dimer scattering in the 
uniform space, f2 = u) = 0. In this case the center-of- 
mass motion separates, /(Ri, R 2 ) = /(Ri — R2), and it 
is convenient to work in the momentum representation, 
/(p) = / /(R) exp(zpR)d 3 i?. Then E c = E c (p) = E — 
p 2 /2/i 3 with fj, 3 = M(M + m)/(2M + to), and 



f(q)/(2n 2 ) d*q 



p 2 + q 2 + (2/i/m)pq - 2pE' 



The calculation of the atom-dimer scattering amplitude 
from Eq. Q then follows the standard method devel- 
oped for the ordinary Schrodinger equation in momen- 
tum space [T9J . In Fig. [I] we present the resulting s- 
and p-wave phase shifts (top) and the corresponding par- 
tial scattering cross-sections (bottom) as functions of the 
atom-dimer collision energy e normalized to the dimer 
binding energy |e | defined by ^/2^|e |— a~ 1 +2fiR*\e \ — 
0. We consider the atom-dimer scattering below the 
dimer breakup threshold, so that the total energy E = 
— |eo| + e < 0. On the left we present the results for 
the mass-balanced case and we see that the s-wave con- 
tribution to the scattering is always dominant. On the 
right, in the case M/m = 6.64, solid lines show the phase 
shifts and scattering cross-sections in the limit of a wide 
resonance, or very small detunings (a 3> R*). One can 
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Figure 1: (color online). The scattering phase shifts and par- 
tial cross-sections as functions of the atom-dimer collision en- 
ergy normalized to the dimer binding energy for the mass- 
balanced case (left) and for the case M/m — 6.64 (right). In 
the top-right graph we also present the dependence on the 
width of the atomic resonance. Dashed, dotted, and dash- 
dotted lines stand for a/R* = 16, 4, and 1, respectively. 
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clearly see that the p-wave contribution dominates ev- 
erywhere except for very small e. The p-wave phase shift 
crosses the unitarity line 6 P = n/ 2 at e ~ 0.1|eo| an d 
stays very close to it at larger collision energies. 

On the top-right in Fig. [I] we present the scattering 
phase shifts in the case of a Fcshbach resonance of fi- 
nite width. We see that the narrower the resonance or 
the larger the detuning, the weaker the atom-dimer in- 
teraction (attractive p-wave as well as repulsive s-wave) . 
This is consistent with the fact that the atoms forming 
the dimer spend more time in the closed channel state, 
interaction of which with the other atom is not resonant. 

At this point we can conclude that for sufficiently small 
detunings R* /a a mixture of K atoms and K-Li dimers 
should behave quite differently compared to the case of 
an atom-dimer mixture in which dimers are composed 
of equal-mass atoms. The resonant character of the p- 
wave atom-dimer interaction can be demonstrated, for 
example, by colliding with each other two cold clouds 
of atoms and molecules and by measuring the angular 
distribution of scattered atoms or molecules [20] . 

We see that the atom-dimer scattering properties are 
sensitive to the detuning R* /a. However, for the Li-K 
mass ratio, increasing this parameter can only weaken 
the atom-dimer interaction. Aside from the scaling with 
| £o 1 1 is it possible to shift the position of the p-wave atom- 
dimer resonance down and eventually turn it into a trimer 
state by changing parameters of the system? According 
to Ref. [15] this can be achieved by increasing the mass 
ratio to the value M/m w 8.2. Another approach, not 
involving changes in atom masses, is based on the fol- 
lowing facts. In the 2D case the K-K-Li system with 
zero-range interactions has one bound trimer state with 
unit angular momentum |21] and in the case of mixed 
dimensions, when K is two-dimensional and Li is three- 
dimensional, the 3-body system exhibits the Efimov ef- 
fect; i.e., it supports an infinite number of trimer states 
[2"2l 125] . Therefore, one can assume that a gradual in- 
crease of an external quasi-2D confinement does shift the 
position of the atom-dimer resonance and turns it into a 
trimer state. 

To answer the question of how strong should be the 
confinement we return to Eq. ([3| and solve it in the 
inhomogeneous case. Note that all the operators in 
Eq. (§ conserve the planar center-of-mass momentum 
as well as the planar angular momentum £ and, there- 
fore, / essentially depends only on three coordinates: 
/(Ri,R 2 ) = f(Z 1 ,Z 2 , |px - p 2 \) exp(i£<p) , where Z l and 
p i are the axial and radial components of R^ , and <f> is the 
angle of p 1 — p 2 . Besides, in the case f2 = uj the center- 
of-mass motion in the axial direction separates and the 
configuration space of Eq. (J3j) becomes two-dimensional. 

Figure [2] shows the trimer formation thresholds for 
uj = f2 (solid line) and uj = (dashed line) , and we intro- 
duce the heavy atom oscillator length, Iq = 1/ y/Mfl. In 
both cases the critical confinement length for R* = is 




Figure 2: Phase diagram of the K-K-Li system in an external 
quasi-2D confinement. The solid (dashed) line corresponds to 
the trimer threshold in the case fi = uj (ui — 0). Three atoms 
form a bound state in vertically (horizontally) shaded areas. 
The dotted line represents the dimer formation threshold in 
the case w = (see text). 

given by lo/a m 1.3 and then rapidly decreases with R*. 
The closeness of the two curves for a < 2l is explained by 
the fact that the dimer size, ~ a, is smaller than the light 
atom oscillator length, loy/MCl/muj, and the dimer and 
trimer states are fairly insensitive to the light atom con- 
finement. Accordingly, this part of the diagram should 
be quite universal for < uj < f2. 

On the negative side of the resonance the two cases 
differ significantly. For finite uj, there always exists a 
dimer state whose size exponentially increases with the 
detuning \lo/a\ (a < 0). The problem then becomes es- 
sentially 2D (cf. [2TJ) and for sufficiently small R* there 
is one trimer state with unit angular momentum (verti- 
cally shaded area). In the case ui = there exists a dimer 
threshold (dotted line) where the dimer size diverges and 
the resulting mixed-dimensional 2D-3D system exhibits 
the Efimov physics [2"2l [25] , We find, however, that the 
corresponding region in the vicinity of this line is 7 or- 
ders of magnitude narrower than the horizontally shaded 
area, which indicates a single trimer state. 

Let us now discuss realistic values of a, R* , and Iq, 
achievable in ongoing experiments. Hereafter we make 
h explicit. For the resonance at B = 114.47(5)G in 
the |l/2,+l/2) Li - 1 9/2, +9/2) K mixture, believed to be 
one of the widest of the 6 Li- 40 K system, the width is 
AB = 1.5(5)G [7J. With the background scattering 
length dbg = 3nm and fi vc i = 1.57 fj,B for the difference 
in magnetic moments of the closed and open channels 
we have [UJ R* = ti 2 /(2^a bg /iroiAB) ~ lOOnm. The 
same value of Iq is achieved for the confinement frequency 

« 2ir x 25kHz. Then the trimer formation threshold is 
reached for a s» 400nm. Note that R e « 2.2nm and the 
zero-range approach is well justified. 

It is well established [24] that atomic systems close to 
a p-wave resonance suffer from very strong losses due to 
the relaxation to deep molecular states. Because of the 
centrifugal barrier the size of a p-wave bound state is of 
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the order of R ei independent of the detuning. Therefore, 
p-wave molecules easily relax to deeper bound states of 
approximately the same size when colliding with another 
atom or molecule. The two-body relaxation rate constant 
is then of the order of a rc i oc hR e /m. Translating this 
picture to our p-wave atom-dimer resonance we observe 
that the range of the atom-dimer potential, given mostly 
by the exchange of the light atom, is one or two orders of 
magnitude larger than R e . However, the relaxation rate 
constant a re i does not scale accordingly as there is only 
a single trimer state and there are no deeper states of 
the same size to which it can relax. Thus, the main loss 
channel in the gas should be the relaxation to molecular 
states of the size R e . This mechanism is responsible for 
the decay of a single trimer, the rate of which we will now 
discuss. It is also important for the atom-dimer mixture 
close to the p-wave resonance as atom triples spend a lot 
of time in this bound or quasi-bound trimer state. 

Consider a K-K-Li trimer in the case R* — and 
a ~ Iq. For hyperradii R e <C £ <C a its wavcfunc- 
tion reads \ff w g"'' -1 /a u * +3 where <&(fi) is a nor- 

malized function of hyperangles and v p « 0.198 corre- 
sponds to the heavy-heavy-light three-fermion problem 
with M/m = 6.64 and p-wave symmetry [TH [THj. The 
decay rate equals 1/r = P(R e )u m i c , where P(R e ) ~ 
(Re/a) 2 ^^ is the probability to find all three atoms 
in the microscopic region £ < R ei and f m i c ~ h/mR 2 
represents the order of magnitude estimate of the re- 
laxation rate once the atoms are confined to this re- 
gion (cf. [23]). In the case of finite R* it follows 
from the coupled-channel analysis that the open-channel 
wavefunction ^(R^Ra,!-) « /(R 1 ,R 2 )/47r|r - R x | is 
always accompanied by the closed-channel amplitude 
•^/i?*/47r/(Ri, R2), and estimates of the normalization 
integral calculated with these functions show that the 
latter is dominant at £ ~ |Ri — R2I <C R* ■ This means 
that at £ <C R* the three-atom problem reduces to the 
problem of a non-interacting atom and a closed-channel 
(bare) molecule. Accordingly, in the case R e -C R* -C a 
the trimer decay rate is a product of three factors, 
1/r = P(R*)v mic (R e / R*) 5 , where the last factor gives 
the probability for a bare molecule and an atom with 
unit total angular momentum to approach each other 
from distances ~ R* to distances ~ R e . Similar anal- 
ysis can be performed for other values of a, R* , and Iq, 
and substituting realistic numbers for Li-K mixtures one 
obtains that the trimer decay rate is 3 to 5 orders of 
magnitude lower than the dimer binding frequency and 
can lead to trimer lifetimes of order Is. In contrast to 
the Efimov case, the binding of these purely long-range 
trimers comes from distances of the order of the dimer 
size, whereas at short distances the three-atom dynamics 
is dominated by the centrifugal barrier. We see that the 
p-wave atom-dimer resonance does not compromise the 
lifetime of the system. 

Independent of the proximity of the p-wave atom-dimer 



resonance the lifetime of the mixture is likely to be lim- 
ited by the relaxation in s-wave atom-dimer and dimer- 
dimer collisions. A detailed discussion of these pro- 
cesses in the case R* = can be found in Ref. [T4"] 
and can be generalized to the case of finite R* in the 
same way as outlined above. In particular, in the case 
R e -C R* -C a the 3D s-wave atom-dimer relaxation rate 
constant equals a rc i » (hR e /m)(R* /a) 5 04 , and we get a 
similar power law dependence for the dimer-dimer chan- 
nel with a slightly lower exponent . We see that for small 
detunings from the Feshbach resonance; i.e., for R* <C a, 
these rates are strongly suppressed. Thus, also in this 
respect the mixture of K atoms and K-Li dimers holds 
promise as a unique long-lived system in which resonant 
p-wave interactions can be observed. 
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